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It is often desirable to determine the potential for radio frequency interference 
between Earth stations and orbiting spacecraft. This information can be used to 
select frequencies for radio systems to avoid interference or it can be used to deter- 
mine if coordination between radio systems is necessary. A model is developed that 
will determine the statistics of interference between Earth stations and elliptical 
orbiting spacecraft. The model uses orbital dynamics, detailed antenna patterns , 
and spectral characteristics to obtain accurate levels of interference at the victim 
receiver. The model is programmed into a computer simulation to obtain long- 
term statistics of interference. Two specific examples are shown to demonstrate the 
model. The first example is a siniu/ation of interference from a fixed-satellite Earth 
station to an orbiting scatterometer receiver. The second example is a simulation 
of interference from Earth-exploration satellites to a deep-space Earth station. 


I. Theory and Models 

Figure 1 contains an illustration of the interference ge- 
ometry for Earth orbiters and an Earth station. Space- 
craft 1 may be transmitting or receiving. Its antenna is 
pointed toward an arbitrary location on Earth. The Earth 
station may be transmitting or receiving. Its antenna is 
pointed toward spacecraft 2 or toward an arbitrary point 
described by the Earth station antenna azimuth and ele- 
vation. 

Two interference scenarios are considered. In the first 
scenario, the Earth station is transmitting a signal toward 
spacecraft 2. This signal is unintentionally received by 
spacecraft 1. In the second scenario, spacecraft 1 is trans- 
mitting a signal that is unintentionally received by the 
Earth station. 

The interference geometry shown in Fig. 1 is common 
to many interference scenarios that occur between two dif- 


ferent radio systems. The level of interference that occurs 
at a victim receiver depends on angles y g and y rg and the 
distance D s j that are shown in Fig. 1. The primary em- 
phasis of the first part of this article is to show how to 
compute these parameters as a. function of time. The an- 
tenna gain of spacecraft 1 in the direction of the Earth 
station is a function of angle y g . The path loss between 
spacecraft 1 and the Earth station is a function of the dis- 
tance D s j. Lastly, the antenna gain of the Earth station 
in the direction of spacecraft 1 is a function of angle y rg . 
The antenna gains and path loss are used to determine 
interference levels at the victim receiver as a function of 
time. 


The interference angles and path distance in Fig. 1 may 
be computed with standard orbit determination methods 
[1]. The appendix contains detailed derivations of the an- 
gles and path distance. 
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The International Telecommunication Union (ITU) an- 
tenna pattern is used to calculate the Earth station an- 
tenna gain [2]. 


G(y rg ) = G P - 2.5 x 1(T 3 



(dBi), 


20x106 (wkn) m <3) 

The angle j g in Fig. 1 is used to compute the antenna 
gain of spacecraft 1 in the direction of the Earth station. 
The spacecraft gain pattern is mission dependent. 


0 deg <~frg <4> m 

~ G i, 4>m ^ yrg < 4*r 

— 32 25 x log (7r^)j 4*r ^ *frg ^ 48 deg 

= — 10, 48 deg < y rg < 180 deg (1) 


where 

G p — peak antenna gain, dBi 
D — antenna diameter, m 
A = wavelength = c//, m 
G\ = 2 + 15 x log (D/A), dBi 
<j, m = 2QX(G P - Gi) 1,2 /D, deg 
<f> r = 15.85(D/A)-°- 6 , deg 
c = speed of light = 3 x 10 s , m/sec 
/ = frequency, Hz 

This equation is valid for D/A > 100. For D/A < 100, a 
different pattern is used. 

G( Jrg ) = G p - 2.5 x 10- 3 (^ 

0 deg < 7 rg < 4 > rn 

_ r , . . „ 100 A 

— G 1 5 t tti ^ 7 rg ^ jj 

= 52-10 log (y) - 25 x log (>„), 



100A 

1) 


< 7 rg < d8 deg 


II. Application 1: Simulation of Interference 
From a Fixed-Satellite Earth Station to a 
Scatterometer Receiver and Calculation 
of Scatterometer Power Flux Density 

A. Introduction 

The World Administrative Radio Conference of 1992 
(WARC-92) allocated the 13.75- to 14-GIIz frequency 
band to fixed-satellite (Earth-to-space) service on a pri- 
mary basis. NASA scatterometers are permitted to use 
this band under the category of space research (spacecraft 
radiolocation) service on a secondary basis. The Interna- 
tional Radio Consultative Committee (CCIR) formulated 
Task Group 7/3 (TG 7/3) to study the use of the 13.4- 
to 14-GIIz frequency band by science services. One of the 
goals of TG 7/3 was to determine the effect that the new 
allocation (to fixed-satellite service) would have on NASA 
scatterometers near 14 GHz. The models developed in 
Section I and the appendix were used to generate a simu- 
lation of interference from a fixed-satellite Earth station to 
a scatterometer in Earth orbit. The scatterometer is rep- 
resented by spacecraft 1 and the fixed-satellite spacecraft 
is represented by spacecraft 2 in Fig. 1. The fixed-satellite 
Earth station is represented by the Earth station in Fig. 1. 

B. Analysis 

The computation of the scatterometer antenna gain in 
the direction of the fixed-satellite Earth station depends 
on components of the angle in Fig. 1 . The scatterometer 
antenna coordinate and spacecraft coordinate systems are 
illustrated in Fig. 2. There are six antenna beams, whose 
positions are defined by angles a, and 7. A vector in 
the antenna coordinate system can be transformed into a 
vector in the spacecraft coordinate system. 


= 10-10 log (y) , 

48 deg < 7 rg <180 deg (2) 

The path loss from spacecraft 1 to the Earth station is 
computed. 


x s = Ax a 



D G 
E H 
F l 


( 4 ) 
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A = cos (7) cos (a) 

B = cos (7) cos (90 — a) 

C = cos (90 — 7) 

D = -cos(90 — a) cos (/?) + cos (90-/?) cos (90 -7) cos (a) 
E = cos (a) cos (/?) + cos (90 — /?) cos (90 — 7) cos (90 — a) 
F = — cos (90 — /? ) cos (7) 

G = -[cos (90— a) cos (90-/?)+cos (9O-7) cos (/?) cos (a)] 
H = cos (a) cos (90 — /?) — cos (90 — 7) cos (/?) cos (90 - a) 
I = cos (/?) cos (7) 

A vector in the spacecraft coordinate system can be 
transformed into a vector in the orbit-plane coordinate 
system. 


where 


= Bx fl 


( 5 ) 


-Mml 

a 

X\v 

a 


L 0 


Uhl 

a 


0 J 


a = semimajor axis (Earth radii) 


fan beam antenna (Fig. 3). The antenna boresight points 
in the k direction. 


cos (/?„) = 


%sf a 


(x>fa 2 + ^/a 2 ) 1/2 


j a 


cos (a„) = , 

(y,fa 2 + Zsfa 2 ) ~ 


( 8 ) 


These angles are used to compute the narrow-beam and 
wide-beam components of the fan beam antenna gain. 


Gl{0g) = -3 (0^73) ’ 0 de S ^ Pt < 0 5 de S 

(narrow-beam component, dBi) 

= -23.9 - 2.16/?^, 0.5 deg < 0 9 < 5.6 deg 
= -48, 5.6 deg < P 9 < 180 deg 

G 2 (a g ) =34-3 (^) 2 , 0 deg < a, < 20 deg 

(wide-beam component, dBi) 

= 35.83 - 0A02a g) 20 deg < a g < 114 deg 
= —14.5, 114 deg < a g < 180 deg (9) 


A vector in the orbit-plane coordinate system can be 
transformed into the right ascension-declination coordi- 
nate system. 

X a n — P x wa (0) 

where P is defined in Eq. (A-3) in the appendix. 

The vector (from the scatterometer spacecraft to the 
fixed-satellite Earth station) is transformed from the right 
ascension-declination coordinate system into the antenna- 
coordinate system. This is done by converting from the 
right ascension-declination system to the orbit-plane sys- 
tem, to the spacecraft system, and to the antenna system. 

x sfa = — A T B P Xfg = {Xsfa j Usja > z sfa) C^) 

These components are used to determine the off-axis an- 
gles in the narrow-beam plane and wide-beam plane of the 


The fan beam antenna gain in the direction of the fixed- 
satellite Earth station is the summation of the narrow- 
beam and wide-beam components. 

G/i = Gi(/?,) + Ga(a,) (10) 

where the minimum value of Gj b is —14 dBi. At this 
point, there is enough information to compute the level of 
interference from the fixed-satellite Earth station to the 
scatterometer. 

/ = P T + G( Jrg) + SL + PL + G Jb (11) 

where 

I = interference power spectral density (dBW/Hz) 
Pt = fixed-satellite Earth station transmitter power 
= 85 - 53.4 dBW 
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G(7r^) — fixed-satellite Earth station antenna gain in di- 
rection of scatterometer [Eq. (1)] 


G p = 53.4 dBi 
D — 4.5 m 
f = 13.995 GHz 

SL = 10 x log (1/BW) (dBW/Hz) 

BW = bandwidth = 2 MHz 

PL = path loss from fixed-satellite Earth station to 
scatterometer [Eq. (3)], dB 


Gfb — fan beam antenna gain in direction of fixed- 
satellite Earth station [Eq. (10)], dBi 


The threshold level where interference is considered 
harmful to the scatterometer is —207 dBW/Hz. 1 The 
ground track of the scatterometer is computed to provide 
a visual indication of where the scatterometer is interfered 
with. The latitude of the ground track is computed first. 


L agt = 90 — cos 1 



( 12 ) 


Then, the longitude is computed by 


360/ 

l ° 9t ~ iol "l436l’ y -° 


■ 


360/ 

1436.1 


, y < o 


(13) 


program computes the maximum and average times of an 
interference episode for the complete simulation period. 
Then histograms of these interference episodes are gener- 
ated and the distribution of interference episodes may be 
estimated. 

The scatterometer measures wind speed over the ocean. 
Taking this measurement requires that three antenna 
beams (beams 1-3 or beams 4-6) have no interference. 
The effect of interference from a fixed-satellite Earth sta- 
tion to the scatterometer is to create certain areas on the 
ocean where measurements of the wind speed cannot be 
taken by the scatterometer. It is possible that the scat- 
terometer will have interference on the ascending portion 
of the orbit over a particular area of the Earth. When the 
scatterometer passes over this area on a descending por- 
tion of the orbit, it is possible that there will be no inter- 
ference. Therefore, it is desirable to compute the statistics 
of interference for areas of the Earth. 

For a simple example, the Earth is divided into areas of 
100 km x 100 km near the equator. To facilitate program- 
ming, this is done by dividing the Earth into areas that 
are 0.89833-deg longitude in width and 0.89833-deg lati- 
tude in height. The width of these areas is 100 km at the 
equator, tapering off to a width of 0 km at the north and 
south poles. The antenna beams are very narrow in one 
direction and wide in the orthogonal direction. Therefore, 
the beams are approximated with 31 vectors in the wide- 
beam plane of the antenna. These vectors are spaced 1 deg 
apart, yielding a total coverage of 30 deg in the wide-beam 
plane. The point of intersection of an antenna vector with 
the Earth is determined. 


where 


j — cos 


X 

a (sin [90 - L agt ]) 


Figure 4 shows the geometry. The unit vector v/|v| 
is converted from antenna coordinates to the right ascen- 
sion declination coordinate system. Then angle a v can be 
calculated. The magnitude of v and the vector R can be 
determined by 


The longitude is transformed to obtain values between 
— 180 and 180 deg. 



a (cos [a v ]) 


C. Statistics and Probability of Interference 

When the scatterometer is visible to the fixed-satellite 
Earth station, it is possible for interference to occur. The 
duration of these interference episodes varies with each 
orbit of the scatterometer. The scatterometer-simulation 


1 M. Spencer, “Interference Criterion for Ku-Band Scatterometer,” 
JPL Interoffice Memorandum (internal document), Jet Propulsion 

Laboratory, Pasadena, California, December 8, 1992. 


|R| 2 = a 2 + |v| 2 — 2a|v| cos (a„) 

R= x + v (14) 

where 

v / \T 

— = (v x ,v yi v z ) 
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This vector location is converted to latitude and longitude 
[Eqs. (12) and (13)]. Then, the area to which this point 
belongs is determined. 

Wind-speed measurements on the right side of the scat- 
terometer are taken by combining signals from beams 1-3 
and, on the left side, by combining signals from beams 4- 
6. An area on the right side is considered to be covered if 
a vector from beams 1-3 intersects the area on the same 
orbit and the interference is below the receiver interfer- 
ence threshold for beams 1-3. Likewise, an area on the 
left side is considered to be covered if a vector from beams 
4-6 intersects the area on the same orbit and the interfer- 
ence is below the receiver interference threshold for beams 
4-6. The probability of interference is estimated as 


The direction of the antenna gain vector is computed 
in the antenna coordinate system (Fig. 3). 

x va = (0,cos [90 - a,], cos [a g ]) T (17) 

This vector is transformed into the right ascension 
declination coordinate system with Eqs. (4-6). 

X V r — ’ P BAx va = (#yr j Vvr ) %vr ) ( 1$) 

Then, the angle (a 5e ) is computed by 

x • x vr = —a [cos (a 5e )] (19) 



(15) The distance between the scatterometer and the point 

of intersection with the Earth can be computed by 


where 

p = estimate of probability of area interference 

n — the number of times that an area is covered and 
interfered with 

d — the number of times that an area is covered 

where the definitions of “covered” and “interfered with 1 
are provided in the paragraph above. 

D. Power Flux Density of the Scatterometer 

Compatibility between terrestrial and spacecraft radio 
systems is often determined from the power flux density 
of the spacecraft on the surface of the Earth. Therefore, 
the power flux density of the scatterometer is computed. 
Figure 5 contains an illustration of the geometry. 


D se = a [cos (a je )] - [(a [cos (a je )]) 2 - (a 2 — 1)] 1/2 (20) 

The vector to the point of intersection of the antenna gain 
vector and the Earth is computed by 

X cr — X + D se X. vr (^1) 

The latitude and longitude of this point, L ae , and L oe) are 
computed with Eqs. (12) and (13). 

It is desirable to know the pow T er flux density as a func- 
tion of elevation angle on the Earth. Therefore, the eleva- 
tion of the vector — x vr is computed. 

Ei = sin -1 (— x vr cos [L ae ] cos [L oe ] - y vr cos [L ae \ 


P ]d = P t + G 2 {c* g ) + FSS (16) x cos [90 - L oe ] - z vr cos [90 - L ae ]) (22) 


where 

Pj d — power flux density of scatterometer at surface 
of the Earth, dBW/m 2 

P t = scatterometer transmit power, dBW 

G 2 (ct g ) = scatterometer transmit antenna gain as a 
function of angle (a^) in the wide-beam plane 
(Fig. 3), dBi 

FSS ~ free space spreading loss, dB/m 2 
= 10 x log (1/[47 t(6378 x 10 3 A, e ) 2 ]) 


E. Simulation 

Table 1 shows the results of a sample simulation. In- 
terference is considered harmful if it exceeds the threshold 
shown on the table. The fixed-satellite Earth station and 
geostationary satellite are arbitrarily located. Interference 
is computed at I-sec intervals. The estimated probabil- 
ity of interference is the number of interference minutes 
divided by simulation minutes. The scatterometer orbital 
period is approximately 100 min. Because the position 
of the orbit changes relative to the fixed-satellite Earth 
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station, some orbits have no interference and other orbits 
have episodes of interference for varying amounts of time. 
The maximum and mean of these interference episodes are 
provided. 


III. Application 2: Simulation of Interference 
from an Earth-Exploration Satellite to a 
Deep-Space Earth Station 

A. Introduction 

The Deep Space Network (DSN) uses the 8400- to 8450- 
MHz band for space-to-Earth transmissions. Earth-station 
receivers are protected by interference criteria that have 
been negotiated in international forums. Other radio ser- 
vices that use this band are aware of the interference cri- 
teria and limit their transmissions accordingly. However, 
radio services that transmit in adjacent bands may not be 
fully aware of their emissions in the DSN band. If these 
out-of-band emissions are strong enough, they can disrupt 
DSN communications. 

B. Power Spectral Density of Adjacent Band 
Interferes 

The power spectral density from an Earth-orbiting 
spacecraft is computed with the following equation: 


Psd — Pt + SD(f) + Gt -f PL\ + Gr (23) 


where 

Psd — power spectral density of interfering space- 
craft at a deep-space Earth station receiver, 
dBW/Hz 


P T — spacecraft transmitter power, dBW 

SD(f) = spectral density of spacecraft transmitter 
(Table 2), dB/Hz 

Gt = peak transmit antenna gain, dBi 
PL\ = path loss, dB 

= 20 x log (c/(47Ti4 m /)] 
c = speed of light = 3 x 10 5 km/s 
A m = minimum orbit altitude, km 
/ = frequency, Hz 

Gr = DSN receive antenna gain = 74 dBi 


The spectral density is computed from the equations 
shown in Table 2 [3]. 2 

In particular, the Earth Observing System (EOS) pro- 
gram is planning the launch of several spacecraft [4] that 
will have spectral emissions in the 8400- to 8450-MHz 
band. Using Eq. (23), the emission of an EOS spacecraft 
is determined to exceed the DSN interference criterion of 
-220.9 dBW/Hz by about 45 dB in the 8400- to 8450-MHz 
band. Equation (23) produces a worst-case power spectral 
density at the deep-space Earth station because it uses 
peak antenna gains and the minimum orbit altitude to 
compute the path loss. Statistics of the interference power 
spectral density are provided next. 

C. Simulation of Interference to a Deep-Space Earth 
Station from Low Earth-Orbiting Spacecraft in an 
Adjacent Band 

It is useful to know the amount of time that the power 
spectral density from an adjacent-band spacecraft exceeds 
the interference criteria of the deep-space Earth station. 
Figure 6 shows a plot of satellite visibility to an Earth sta- 
tion versus orbit altitude. It is assumed that the satellite 
passes directly over the Earth station with a circular or- 
bit. The geometry is illustrated in Fig. 7. Visibility time 
is computed as 


^ _ 36001iff,_ N 

Pvis — (sec) 


(24) 


where 

T 0 — orbital period, hr [5] 

= 2tt(^ + ^) 3/2 /m 1/2 

A = orbit altitude, km 
R = Earth radius — 6378 km 
fj. = Earth gravitation, mass product 
= 5.17 x 10 12 km 3 /hr 2 
9 = central angle, rad 
= cos _1 (R/[R -j- A]) 


If the satellite has sufficient power, it can produce line- 
of-sight interference to the Earth station during its orbit 
visibility time. This can be a significant amount of time 
even for a low-altitude spacecraft. Interference times for 


2 B. Tunstall, “Transmit Filtering for EOS Direct Access System 
Downlink,” Memorandum to J. Scott (internal document), God- 
dard Space Flight Center, Greenbelt, Maryland, July 10, 1992. 
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actual spacecraft systems depend on transmitter power, 
orbits, antenna gains, etc. Figure 1 can be used to il- 
lustrate the interference geometry. The EOS spacecraft 
is represented with spacecraft 1. The deep-space Earth 
station points at a specified azimuth and elevation. The 
EOS spacecraft is in a circular orbit around the Earth with 
its antenna pointing toward the center of the Earth. The 
angle j rg between the boresight of the deep-space Earth 
station and the vector to the EOS spacecraft is computed 
with Eqs. (A-15) and (A-17) in the appendix. The deep- 
space Earth station (70-m) antenna gain in the direction 
of the EOS spacecraft is computed with Eq. (1). 

G(j rg ) = 74 - 0.0025 (1960 7r5 ) 2 (dBi), 

0 deg < 7 rg < 0.0485 deg 
= 51.4, 0.0485 deg < j Tg < 0.168 deg 
= 32 - 25 log (7^), 0.168 deg < j rg < 48 deg 
= - 10, 48 deg < j r 9 < 180 deg (25) 

The angle j g between the antenna boresight of the EOS 
spacecraft and the vector to the deep-space Earth station 
is computed with Eqs. (A-19) and (A-20) in the appendix. 
The EOS spacecraft antenna gain in the direction of the 
deep-space Earth station is computed. The model is a 
set of straight line segments that approximates the actual 
pattern. 3 

G e { y g ) = —3.5+ ^y^(dBi), 0 deg < 7s < 3 deg 

= 0.091 - 3 deg < y g < 14 deg 

14t 

— ~ 8 - 5 + 14 deg < 7g < 70 deg 

5o 

= 84.25 - 70 deg < j g < 90 deg 

= -14.5, 90 deg < j g < 180 deg (26) 

The interference power spectral density level at the deep- 
space Earth station receiver is computed by 

I 0 = SD(f ) + P e + G,{ 7j ) + PL + G(y rg ) (27) 

3 J. Hart, “Interference Analysis for EOS AM Direct Access System,” 
Interoffice Memorandum (internal document), Stanford Telecom, 
Seabrook, Maryland, July 7, 1992. 


where P e = EOS transmitter power =11.8 dBW. 

The first simulation on Table 3 is for one EOS space- 
craft transmitting in the direct playback mode (QPSK, 
150 MSPS per channel). Statistics of the interference are 
provided. The second simulation in Table 3 is for five dif- 
ferent spacecraft. The antenna pattern for the last four 
spacecraft is the same as for the EOS antenna [Eq. (26)] 
except for additive constants to provide the correct peak 
antenna gain for each spacecraft. Each interference event 
is composed of a number of interference samples that oc- 
cur at 0.5-sec intervals. Samples are plotted on Fig. 8 for 
the second simulation in Table 3. At each sample interval, 
the amount of interference exceeding the —220.9 dBW/Hz 
criterion [6] is plotted. Some events have only one sample 
(0.5 sec of interference) and some have up to 20 samples 
(10 sec of interference). The simulation ran on a Sun work- 
station. 

The antenna pattern of the deep-space Earth station 
provided in Eq. (25) was developed a number of years 
ago for generalized Earth stations. More accurate gain 
data have been made available for the 70-m antennas of 
the DSN. 4 These data have been fitted with a number of 
equations. 

Gmilrg) = 74.15 - 0.0025(2400 7rp ) 2 (dBi), 

0 deg < -) rg < 0.0376 deg (dBi) 

= 53.7, 

0.0376 deg < y rg < 0.04 deg 
= 57.4 - 0.025(1350[ 7t .£ - 0.049]) 2 , 

0.04 deg < j rg < 0.0626 deg 
= 49, 

0.0626 deg < y rg < 0.0905 deg 
= 25-23 log ( 7r ,), 

0.0905 deg < y rg < 33.2 deg 

= - 10 , 

33.2 deg < 'frg < 180 deg (28) 

4 D. Bathker, personal communication, DSN International Fre- 
quency Management Program, TDA Planning Section, Jet Propul- 
sion Laboratory, Pasadena, California, October 1993. 
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Figure 9 contains plots of antenna data, curve fits from 
Eq. (28), and the ITU pattern [2]. Equation (28) fits the 
data quite well for this application. The ITU model overes- 
timates the gain data in every region except that from 0.04 
to 0.06 deg. It is expected that the interference from low- 
Earth-orbit spacecraft would be reduced when the DSN 
Earth station antenna is modeled with Eq. (28). 

The computer simulation was conducted with G m (j rg ) 
replacing G(y r g ) m Eq. (27). In the first simulation, one 
EOS spacecraft transmitted in the direct playback mode. 
In the second simulation, five different spacecraft were 
transmitting. The interference percentage is one-fourth of 
that in Table 3 with the new antenna model for one EOS 
spacecraft. For five spacecraft, the interference percent- 
age is approximately one-fourth of that in Table 3 with 
the new antenna model. 

All of the simulations use an elevation of 5 deg and 
an azimuth of 100 deg for the deep-space Earth station 
antenna. When larger elevation angles are used, the in- 
terference percentages are decreased. This is because the 
spacecraft spends less time within the same Earth station 
antenna beamwidth at higher elevations. The interference 
levels are computed at 0.5-sec intervals in all simulations. 

IV. Summary and Conclusions 

A model is developed that will determine the statistics 
of interference between Earth stations and elliptical orbit- 
ing spacecraft. Detailed models of the geometry of Earth 
orbiters and Earth stations are shown. It is possible to 
determine antenna gains and path distances quite accu- 
rately with these models. Interference levels at the victim 
receiver are determined by antenna gains, path distances, 
and radio system characteristics. The models may be used 
to predict interference in a wide variety of scenarios that 
involves Earth stations and Earth orbiters. Two examples 
that make use of these models are shown. 


The first example is a simulation of interference from 
a fixed-satellite Earth station to an orbiting scatterome- 
ter receiver. Table 1 contains a summary of a simulation 
that was performed. Models were developed to estimate 
the probability that certain areas of the Earth would not 
have accurate measurements due to interference. Also, 
the power flux density on Earth of the scatterometer sig- 
nal as a function of elevation angle is determined. Results 
from the simulation showed that the scatterometer should 
use a different frequency to avoid interference from several 
fixed-satellite Earth stations. Simulation results compared 
favorably with other simulation programs. 5 

The second example is a simulation of interference from 
an Earth-exploration satellite to a deep-space Earth sta- 
tion. Worst-case levels of power spectral density from low 
Earth orbiters at the station are computed. These lev- 
els exceed the interference criteria of the deep-space Earth 
station in the 8400- to 8450-MHz band. A simulation of 
interference from low Earth orbiters to a deep-space Earth 
station is conducted. This simulation computes the path 
loss and off-axis antenna gains as a function of orbital po- 
sition of the low Earth orbiter. It can be used to predict 
the statistics of interference to the deep-space Earth sta- 
tion. Table 3 contains a summary of two different sim- 
ulations that were performed. Results from the simula- 
tion showed that excessive coordination could be avoided 
if Earth-exploration satellites reduced their emissions by 
20 dB in the 8400- to 8450-MHz band. The interference 
statistics are in close agreement with other simulations. 6 

The simulations are programmed in Fortran and run 
on a Sun workstation. The author may be contacted for 
instructions on the operation of the simulation programs. 


5 J. Hart, personal communication, Stanford Telecom, Seabrook, 
Maryland, April 1993. 

6 J. Hart, “Compliance of EOS Direct Access System with DSN Pro- 
tection Criteria,” Interoffice Memorandum (internal document), 
Stanford Telecom, Seabrook, Maryland, May 17, 1993. 
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Table 1. Example of an Interference simulation from a fixed-satellite Earth station 
to an orbiting scatterometer receiver. 


Parameter 


Value 


Scatterometer interference threshold, dBW/Hz 
Scatterometer altitude, km 
Fixed-satellite Earth station 


Geostationary satellite 
Simulation minutes 

Number of interference computations per minute 
Antenna beam element 
Number of interference minutes 
Estimated probability of interference 
Maximum period of an interference episode, nun 
Mean period of an interference episode, min 


-207 

797 

85-dBW EIRP, 2-MHz bandwidth 
4.5-m antenna 
34 deg N, 128 deg W 

196 deg W 
10,000 
60 
1 

2.067 

2.067 x 10" 4 
0.567 
0.129 


Table 2. Spectral density equations. 


Modulation Spectral density (baseband)* 

Quadriphase-shift keying [2/ SR) X [sin (2 tt / / S R)] 2 / (2tt / / S R) 2 

Unbalanced quadriphase-shift keying (77/Sftj) X [sin (7r / / SRj)] 2 /(tt / / SR/) 2 

+ ( r Q /SR Q ) X [sin (*f/SR Q )] 2 /(7rf/SR Q ) 2 

Phase-shift keying (1/SR) X [sin (tt/ /SR)] 2 /(tt / /SR) 2 

Minimum-shift keying [16 /(tt 2 SR)] X [cos (2tt f / S R)] 2 / (1 — 1 6/ 2 /5/? 2 ) 2 


SR = 
/ = 
ri = 
SR / = 

r Q ~ 

SRq = 


total output symbol rate 
frequency* 

ratio of power in I channel to total power 
symbol rate of I channel 
ratio of power in Q channel to total power 
symbol rate of Q channel 


Replace / with / - f c (where / c is the center frequency) to obtain the spectral 
density at the center frequency. 


Table 3. Simulation results comparing one EOS spacecraft to five other 
spacecraft— EOS (DP), SPOT-4, IRS-1 B, POEM 1, and RADARSAT-1. 


Parameter 

Duration of simulation, min 
Interference events 
Interference percentage 
Interference duration, sec 
Shortest 
Longest 
Average 

Time between interference events, days 
Shortest 
Longest 
Average 


EOS spacecraft 

Five spacecraft 

525,600 (1 year) 

525,600 (1 year) 

37 

128 

8.89 X 10“ 4 

2.39 X 10~ 3 

3.5 

0.5 

10 

10 

7.58 

5.90 

1.99 

0.08 

13.44 

12.83 

10.11 

2.87 



Fig. 1. Interference geometry for Earth orbiters and an Earth 
station. 




Fig. 3. Off-axis angies in the wide-beam and narrow-beam planes. 



Fig. 4. intersection of antenna beam 
with Earth. 


SCATTEROMETER 



Fig. 5. Geometry for scatterometer power flux density 
calculation. 



101 102 10 3 10 4 105 

ORBIT ALTITUDE, km 

Fig. 6. Satellite visibility to an Earth station versus orbit altitude 
for circular satellite orbits. 
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Appendix 

Interference Geometry and Orbital Dynamics Models 


The spacecraft orbit plane is illustrated in Fig. A-l. 
The position of spacecraft 1 in its orbit plane is computed 
first. 



Q x 

Qy 

Qz 


w x ~ 

W* 


x w = (a[cos ( E )] — ae, a[l — e 2 ] ly ^ 2 [sin (E)] t 0) T (A-l) 

where 

a = semimajor axis, Earth radii 
E = eccentric anomaly, rad 
M — E — e(sin [E]) = mean anomaly, rad 
= nt 

e = orbit eccentricity 
n — 0.07437 /a 3 / 2 , rad/min 
t = time, min 

Newton’s iteration [7] may be used to solve for the ec- 
centric anomaly, as follows: 


P x — cos (w) cos (Q) — sin (w) sin (Q) cos (i) 

P y = cos (cj)sin (fi) + sin (u>) cos (fl) cos (t) 

P z — sin (w)sin (?) 

Q x — —sin (w)cos (Q) — cos (uz)sin (Q)cos (z) 
Q y — — sin (u;) sin (Q) 4- cos (w) cos (Q) cos (z) 
Q z = cos (w) sin (*) 

W x — sin (Q) sin (z) 

W y = — cos (fi) sin (z) 

W z = cos (0 
u = argument of perigee 
Q = longitude of the ascending node 
i = orbital inclination 


e(sin [M]) 

E ° M i _ s i n -f e ] -j- sin [M] 
M k = E k - e(sin [E k ]) 


The location of the Earth station is determined in the 
right ascension-declination coordinate system. 

x g — (sin [0] cos [0],sin [0] sin [0],cos [0]) T 


Ek + i = 


Ek + 


M -M k 
1 - e[cos (E k )] 


(A-2) 


where 


k - 0 , 1 , 2 ,... 

The position of spacecraft 1 in the orbit plane is converted 
to its position in the right ascension-declination coordi- 
nate system. 


X = Px w - (x, y, z) T (A-3) 


where 


( x g i Vg j z g ) 


(A-4) 


where 

0 = 90 - l a 

6 — / o + fc + 360*/1436.1 

k — an arbitrary constant used to rotate the Earth 
station relative to the orbit plane 

l a = Earth station latitude 

l 0 = Earth station longitude 

The Earth station may be pointed at another satellite 
(spacecraft 2) or its boresight direction may be described 
with azimuth and elevation notation. If it is pointed at 
spacecraft 2, that satellite’s position may be described 
with Eqs. (A-l) to (A-3). 
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X2 = (z2,y2,Z2) T (A-5) 

A special case exists when spacecraft 2 is in a geosta- 
tionary orbit. It may be assumed that it is located in the 
equatorial plane. The location is determined in the right 
ascension-declination coordinate system. 

*2 = x gg = 6.6257(sin [0^] cos [6 g \ ) sin [0^]sin [0 5 ],O) T 

(A-6) 


= cos 1 (^^ (deg) 

" ,= cos "' (^7) +t+ THir- s '“ s 0 < d '®) 

= 36 °- cos ' , fe ) + t + mi - ! '" <0 

(A-8) 


where 

0^ = 90 

0g = geostationary longitude +k -f 360^/1436. 1 

If the Earth station is pointed at a specific elevation and 
azimuth, its boresight direction may be converted from 
azimuth-elevation coordinates to Earth-centered coordi- 
nates as 

^ec — 1 *— (%ec j Vec j z ec) (A-7) 

where 


where 

r ec = (*ec + Vec + *«) 1/3 
r xy = (*ec + Vic ) 1 ' 2 

Finally, these coordinates are converted to the right 
ascension-declination coordinate system by 

^bor — Tec (sin [0*] cos sin [0ft] sin [0ft], cos [0 6 ] ) T 

(A-9) 


Cl 

c 2 

C 3 

C4 

c 5 

Ce 

c 7 

C 8 

C 9 


Ci — cos (/ 0 ) cos (90 — l a ) 

C 2 — — cos (90 - / 0 ) 

C 3 = cos (la) cos ( l 0 ) 

Ca = cos (90 — l 0 ) cos (90 - / a ) 

C5 = cos ( l 0 ) 

Ce = cos (/ a ) cos (90 - l 0 ) 

Cj — — cos (/ a ) 

c 8 = 0 

C9 = cos (90 - l a ) 

x ae = (—cos [el] cos [az],cos [e/]sin [az],sin [e/]) T 
el = elevation of Earth station antenna 
az = azimuth of Earth station antenna 

These coordinates are converted to spherical coordi- 
nates, as follows: 


The line-of-sight visibility of the Earth station to space- 
craft 1 is determined. Figure A-2 illustrates the central 
angle between the two. The central angle between space- 
craft 1 and the Earth station is computed as 

x • x g = |x| cos (7) (A-10) 

where 7 = central angle between spacecraft 1 and the 
Earth station. 

Figure A-2 also illustrates the limit of visibility. The 
constraint is that the line between spacecraft 1 and the 
Earth station is tangent to the Earth. Central angles that 
are less than or equal to this angle indicate that spacecraft 
1 is visible to the Earth station. 

7 - = c “"'(ki) <A - n) 

The visibility condition is stated as 

7 <7, (A- 12) 
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Interference between spacecraft 1 and the Earth station 
can occur only if they are visible. If spacecraft 1 is visible, 
additional computations are necessary. The angle j rg on 
Fig. 1 is used to compute the antenna gain of the Earth 
station in the direction of spacecraft 1. The vector from 
the Earth station to spacecraft 2 is computed as 

Xf g = X 2 — Xg (A- 13) 

A special case occurs when spacecraft 2 is a geostation- 
ary satellite, as follows: 

x 2 = x gg (A- 14) 

If the pointing of the Earth station antenna is described 
in terms of azimuth and elevation, then 

Xfg = Xbor (A- 15) 

The vector from the Earth station to spacecraft 1 is 
computed as 

Xf s = x — Xg (A- 16) 


The angle between these two vectors [Eqs. (A- 13) and 
(A- 16)] can be used to determine the antenna gain of the 
Earth station in the direction of spacecraft 1. 

Xfg •Xfs = | Xfg j |Xf s | COS (j rg ) (A-17) 

The distance between the Earth station and spacecraft 
1 is computed as 

D s j — 6378|x g — x|(fcm) (A-18) 

The angle between the antenna boresight of spacecraft 
1 and the direction to the Earth station is computed as 

Xbi*(-X fs ) = |x b i||x fs |cos ( 7 g) (A-19) 

where Xbi = antenna boresight vector of spacecraft 1 . 

A special case occurs when the antenna of spacecraft 1 
points toward the center of the Earth. 

x b i = -x (A-20) 
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